Transmission of correlated electrons through sharp domain walls in magnetic 
nanowires: a renormalization group approach 
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The transmission of correlated electrons through a domain wall in a ferromagnetic one dimensional 
system is studied theoretically in the limit of a domain wall width smaller or comparable to the elec- 
tron Fermi wavelength. The domain wall gives rise to both potential and spin dependent scattering 
of the charge carriers. Using a poor man's renormalization group approach for the electron-electron 
interactions, we obtain the low temperature behavior of the reflection and transmission coefficients. 
The results show that the low-temperature conductance is governed by the electron correlations, 
which may suppress charge transport without suppressing spin current. The results may account 
for a huge magnetoresistance associated with a domain wall in ballistic nanocontacs. 

PACS numbers: 73.63.Nm, 71.10.Pm, 75.70.Cn, 75.75.+a 



I. INTRODUCTION 

Domain walls (DWs), i.e. the boundaries separating 
different domains of homogeneous magnetizations^ are 
recently a subject of extensive theoretical and experimen- 
tal investigations. This renewed interest in DWs is stim- 
ulated by their possible applications in magnetic logic 
elements and other nanoelectronics and spintronics de- 
vices. Two effects associated with DWs are of particular 
interest. The first one is the way a DW affects electronic 
transport, i.e., the associated magnetoresistance. The 
crucial point here is that the influence of a single DW on 
the resistance can be controlled by an external magnetic 
field 2 * 3 -. The second effect concerns the influence of elec- 
tric current on the DW behavior (DW motion, magnetic 
switching)!, which allows controlling of DWs by means 
of an electric field. 

Recent advances in experimental techniques have made 
possible the determination of the resistance of a sin- 
gle DW in submicron structured samples^^iLSiSiiSiii. 
The results on the single domain resistance are differ- 
ent in magnitude and sometimes differ also in the sign. 
In the case where the DW width is large on the scale 
set by the Fermi wavelength, Xp, of the carriers, the 
theory of the DW contribution to electrical resistance 
is well-established 12 i 13 i 14 i 15 i 16 i 17 . The spin of the elec- 
tron moving across the wall changes its orientation quasi- 
adiabatically (or even adiabatically for very thick DWs) . 
However, the DWs formed at nanoconstrictions may be 
atomically sharpA^iiSiSSi^i and the spin of an electron 
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crossing the wall does not change quasi-adiabatically. Ac- 
cordingly, a completely different approach to the trans- 
port theory through DWs is required. This is particularly 
true for ferromagnetic semiconductors which are consid- 
ered to be most promising for spintronic applications 22 . 
Indeed, recent experiments on magnetic nanostructures 
and nanowires indicate that the presence of DWs may 
result in a magnetoresistance (MR) as large as several 
hundreds&S 3 . or even thousands 2 ^ 5 - of percents, as op- 
posed to the case of thick on scale of (or adia- 
batic) DWs in bulk metallic ferromagnets. In the bal- 
listic regime, the theoretical treatments towards explain- 
ing this effect 26 i 27 i 28 i 29 i 30 ' 31 rely on the assumption that 
the DW is sharp enough to be treated as a spin depen- 
dent scatterer for the charge carriers. The success of 
these theories in explaining the extraordinary large MR 
is moderate, in particular for metallic ferromagnets such 
as Ni where some features of the physics governing the 
behaviour of MR are still unclear—. 

Another feature of the DWs created at nanoconstric- 
tions is their small lateral size (cross-section of the con- 
striction). This small size limits the number of quantum 
channels active in transport to a few ones or even to 
a single one. Consequently, the constriction behaves as 
a one- or quasi-one-dimensional system. In such a case, 
the role of electron-electron interactions may be crucial 33 
for understanding the behavior and basic transport char- 
acteristics of the DWs formed at nanoconstrictions. It 
is well-established that electronic correlations in a one- 
channel wire result in a non-Fermi-liquid behavior - thus 
forming a Luttinger liquid^a^ 5 -. It is also known that an 
impurity present in the ID Luttinger liquid suppresses 
the linear conductivity, which vanishes even for a weak 
impurity scattering potential 3 ^ 7 -. This can be traced 
back to a vanishing density of states at the Fermi level. 
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At finite applied voltages the transport through the wire 
does not vanish due to the nonlinearity of the current- 
voltage characteristics^. Since a sharp domain wall acts 
in a one-channel wire as a localized spin-dependent scat- 
tering center, one can expect a strong influence of elec- 
tron correlations on the MR at low temperatures. 

To confirm this theoretically one could use bosoniza- 
tion techniques". However, we will follow an- 
other route based on the "poor man's" renormalization 
method^S^i^. In our case, the DW scatters both the 
charge and spin of the carriers. As shown below, our 
scheme allows us to obtain results for the renormalized 
transmission and reflection coefficients in terms of the 
uncorrelated spin-dependent ones (i.e., in terms of the 
reflection and transmission coefficients of the wall in the 
absence of electron-electron interactions). The uncorre- 
lated quantities can be obtained from other schemes, such 
as the Hartree-Fock or density-functional theory (within 
local density approximation) and then used as an input in 
our results to obtain renormalized transmission through 
the DW. Hence, our approach - in combination with nu- 
merical (effective single particle) methods - offers a new 
possibility to understand the material-dependent MR as- 
sociated with a DW creation (destruction) , and possibly 
to resolve some controversy concerning huge magnetore- 
sistance in some ballistic nanocontacts. 

The paper is organized as follows. In Section^ we in- 
troduce the problem and the non-interacting scattering 
states for a sharp domain wall. In Section II I II we use 
perturbation theory in the electron-electron interaction 
to calculate corrections to the scattering amplitudes. We 
obtain the renormalization group differential equations 
for the scattering amplitudes. In Section Hvl we describe 
the zero temperature fixed points predicted by the scal- 
ing equations and the power law behavior of the reflec- 
tion and transmission coefficients of the DW as T — > 0. 
In Section E] we discuss the relevance of our findings to 
realistic physical systems and summarize our results. 

II. MODEL 

We consider a magnetized system with electrons be- 
ing constrained to move in one dimension while being 
exchange-coupled locally to the space varying magnetiza- 
tion, M(r). The wire itself defines the easy (z) axis, and 
a domain wall centered at z = separates two regions 
with opposite magnetizations, M z (z — ► ±oo) = ±Mrj. 
Assuming M(r) to lie in the xz plane, and the domain 
wall to be thinner than the Fermi wavelength, we write 
the single-particle Hamiltonian as: 

Ho = ~!i-^2+ hVS ( z ) + J M z(^z + h\5(z)a x , (1) 
2m dz 

where the term h\6(z)a x describes spin scattering pro- 
duced by the M x (z) component^ 

A = — / M x {z)dz, 

J—oo 



and V is a potential scattering term. Single electron 
wavefunctions are spinors with components Xo{ z ) satis- 
fying the condition 

" 5T- ( J-X*(0 + ) - !-*r(0-)) +V X <r(V)-*X-<r(0) = • 

2m \ oz az I 

(2) 

The electron's wavevector in each domain is related to 
the energy E by 

k=^(E±JM ). (3) 

The electron gas in the negative semi-axis (z < 0) is pre- 
dominantly j-spin. An electron incident from the left 
with the momentum k and spin f (or j) can be transmit- 
ted to the positive semi-axis while preserving its spin, 
but the energy conservation requires the momentum to 
change from k to fc_ (or k+) defined by: 



k ± = Jk 2 ±^ T JM . (4) 

If the transmission occurs with spin reversal, the momen- 
tum k is not changed. 

domain 
wall 

Wire: \ | / ' +~ ~ 

I 

x 



Spin parabolic bands: 




FIG. 1: Schematic of a domain wall and the relevant electron 
spin bands. States lp p f and VV-I have the same energy. 

We label the states through the incident wave, so that 

Mz)= y r '(k)e-* k - z ) ' Z< ° (5) 

describes a scattering state with a wave incident from 
z = — oo with spin f and momentum k > 0. Reflection 
amplitudes of a spin a electron with or without spin re- 
versal are denoted by r' a and r a , respectively. The same 
convention applies to the transmission amplitudes t' a , t a . 
The transmitted wave corresponding to Eq. (JSJ is 

^)=(^2)^)' Z> ° ^ 
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and the scattering amplitudes are given by: 
2(v + u_ +2iV)v 



f T (fc) 
t\(k) 



(v + v_ + 2iV) 2 + 4A 2 
AiXv 

(v + v- +2iV) 2 + 4:\ 2 



r t (fc) + l, (7) 



= r> t (*), 



(8) 



where we have defined the velocities vr±\ — hkr±)/m. 

The scattering state corresponding to a wave incident 
from the left with |-spin is: 



V>m( z < o) = 
i) Ki {z > o) = 



ik . z 



\{k)e 
e"°" + rj_(fc)e 

t\{k)e lkz 
h(k)e ik + z 



and the corresponding amplitudes are: 
2{v + v + + 2iV)v 



t[(k) 



(v + v + + 2iV) 2 + 4A 2 
4i\v 

(v + v+ + 2iV) 2 + 4A 2 



ri(k) + 1. 
rl(fc). 



(9) 

(10) 
(11) 



The expressions for the scattering states corresponding 
to the waves incident from +oo are: 



and 



ip-k,t(z < 0) 
i/i-k,t(z > 0) 

i/>-k,l(z < 0) 
i/)-k,i(z > 0) 



H(k)e "' 



— 

—ikz 



-, — ikz 



r[(k)e lk + z 



t\{k)e~ ikz 

r\{k)e lk - z 
e~ lkz +r ] {k)e lkz 



(12) 



(13) 



where we consider k > 0. We shall henceforth denote by 
e(=tp, |) the eigenenergy of a scattering state with mo- 
mentum +p (or — p) incident from the left (or right). 
The scattering amplitudes satisfy some general relations 
that can be found from a generalization of the Wronskian 
theoremM to spinor wavefunctions. We provide such re- 
lations in Appendix 1X1 

In order to deal with the electron interactions, it is con- 
venient to rewrite the scattering states in second quan- 
tized form, making use of right (d qa ) and left (b qa ) mov- 
ing plane- wave states. 

The operators for the scattering states with electrons 
incident from the left (cfe jCr ) are: 



dq 
2^ 



+ 



q — k — iO q — k— a + «0 



<(k) t if*(k) . ir'*{k) - 

q + k~l0 ga+ q-k + i0 aq '- a q + k_ a -lQ ^ 

(14)" 



and the operators for scattering states with electrons in- 
cident from the right (dk >a ) are: 



dk,a = I — 



2tt 



q + k + iO q + k a — iO 



ir*_ a (k) 



it'*Ak) 



ir'* a (k) 



' q-k + lO qa q + k~i0 ' q-k a +i0 ^ 

(15) 

where denotes a positive infinitesimal and the k- 
subscript a = ±1. By inverting these equations, we ob- 
tain the plane wave operators as linear combinations of 
the scattering-state operators: 



dk 
2ttI 



Ck,, 



t-a{k)c k+ , 



k — p — iO fc — p + iO 



t'_ a (k)c k -q r- a (k)dk,<T _ r'_ <T {k)d k+ -i 
k — p + iO k — p + iO fc — p + iO 



(16) 



dk 
2nd 



dk,a t a {k)d k _ i<y 
k — p — iO k — p + iO 



t' a {k)d k - a r a (k)c k ^ _ r' a {k)c k _ 
k — p + iO k — p + iO 



P 



iO 



(17) 



The electron interactions are introduced in the Hamil- 
tonian through the term: 



Hi, 



dkidq 



at It 



int - 9l,a,0 J (2tt)2 a k 1 , a k2^ ak 2+q,l3 ki-q,a 

/dk\dq * ft f 
(2?r) 2 a kl,a b k 2 ,/3 b k2+q,l3 a k 1 -q, a .(18) 

The coupling constants g\ and 172 describe back and for- 
ward scattering processes between opposite moving elec- 
trons, respectively. The Greek letters denote here the 
spin indices, and the summation over repeated indices 
is implied. Because the Fermi momentum depends on 
spin, we allow for the dependence of g on the spins of the 
interacting particles. We therefore distinguish between 
9it > 9il 9i± and 92h 92[, 92±- The forward scattering 
process between particles which move in the same direc- 
tion will not affect the transmission amplitudes, although 
it will renormalize the Fermi velocity^. This effect is 
equivalent to an effective mass renormalization and the 
electrons with different spin orientations may turn out to 
have different effective masses, in which case our calcu- 
lations remain valid, as shown in Appendix 151 



III. SCALING EQUATIONS 

The corrections to the transmission amplitudes will be 
calculated to first order in the perturbation Hi nt . It has 
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been shown in Ref. 40] that the corrections diverge log- 
arithmically near the Fermi level. These divergences will 
later be dealt with in a poor man's renormalization pro- 
cedure. 

Let us consider the Matsubara propagator, 



dQ H(-Q X)r\{Q) _ f(Q T)^(Q) ' 

2ni V 2Q - q - iO 2Q-q + iO 



(22) 



where f(iQ T) denotes the Fermi occupation number of 
the state tp±Q^. In order to calculate the propagator 
{T T a Pi |(r)6{, 2 ^) we again expand ^ using Eq. 117L 



-ffi (t')cLt' t ' K3,\i o "2T "fc 2 ,T 

0(r) = — (T T e J T %>,t( t )v,-|-}o : (19) and then with help of Eg. HA7fl we obtain 



where (...)o denotes the average in the non-interacting 
Fermi sea. The propagator for non-interacting electrons 
is then given by: 



S (0) M = 



e(p', T) 



p — p' + iO p — p_ — iO 



The transmission amplitude appears associated with the 
denominator p — p'_ — iO which, for the variable p, gives 
a pole in the upper half plane. The meaning of this pole 
is that the transmitted particle is a right-mover in the 
z > half-axis. Our strategy is to calculate the first 
order correction term (in Hi nt ) to Q, which will have 
the same form as the second term in (|20() . so that the 
amplitude correction, 8t^(p'), can be read off from the 
result. We now explain the procedure in some detail. 

We begin by considering the first order expansion for 
Q in the coupling g\^. For simplicity, we shall henceforth 
omit the subscript "0" in the brackets, since we will be 
dealing with the non-interacting Fermi sea, unless other- 
wise stated. From Wick's theorem we get the first order 
correction to the propagator in equation l|20l) as 

S (1) (t)=3i T [ 
(a P ,T ( r )&l 2 , T ( T ')) (4 lt T ( T ')^ fe i -9.T ( r ')> (ofc 2 +9,T ( T ')<v, T ) 

+ (apj(r)oj CiT (T0)(6i 2iT (T / )afe 2 + g ,T(T / ))(fcfe 1 - g ,T(r')c| ) ',T> 

(21)' 

where the internal momenta kin), q and time r' are to 
be integrated over and the time ordering T T is implicit. 
There are also two other Wick paired terms at instant r' 
of the form (a* (t>(t')} and (ft (T J )b(T J )). We have omit- 
ted these terms in Eq. I|21l) because they will not be loga- 
rithmically divergent: the divergences arise from electron 
reflection by the Friedel oscillations in the Fermi sea^. 
Such reflection processes appear in equation (|21|l through 
(St(r')a(r')> and (at( r ')S( r ')). 



To calculate 



'3 f . 

*i,T 



we make use of the ex- 



pression i|17|) for bk-L-q,]- The contour integration over 
k\ eliminates the terms containing poles in the same 
half-plane. Fermi sea averages, such as (a\. x -t-c^-f) and 

(a\. x fdkrf), can be calculated in the same way as in 
Eq. 'fZOjL The result is: 



dki t - 



l/T 



dQ 
2tt 



dre 1 



1 



1 



- e(— Q T) 



p' + iO p-Q-iO 



1 



rUQ) 



iuj - 



e(Qt) 



p' - iO p-Q + iO 



(23) 



The presence of two different energy poles can be under- 
stood from the fact that a qa (or b qrT ) represents a plane 
wave running over the entire z axis and its energy cannot 
be the same on both sides of the domain wall because of 
the energy dependence on spin. 

Using Eqs. 12 I'll and 12MI) we can calculate the first term 
in Eq. ||3T} as 

(apj(T)^ 2:T (T'))(4 1 , T ( r ')^fei-g,T(' r ')>(a/c 2 +g,T(T / )Cp/ iT ) 



dQidQ 2 



e{p'])J (2^) 2 iuj-e(-Q 2 ]) 
ri(Q2)rl(Qi)t* r (p')f(-Qt T) 



(P 



»0)(2Qi 



iO) 



(24) 



The analytic continuation of the Green's function fre- 
quency, iuj — > ui + iO, gives the retarded Green's function. 
The frequency denominator (iu> — e(~Q 2 t)) yields a 
principal Cauchy part plus a delta function part. The 
latter isolates the energy pole at e(— Q 2 T) — w and we 
choose u> = e(p' f) => Qi = v'-- We shall only retain this 
delta function part. Therefore, we put Q 2 = p'_ in the 
integrand and, by comparing with (|20|l , we conclude that 
the contribution of the first perturbative term in Eq. I|21|) 
to the transmission amplitude is given by 



SiT f dQ f(-Q t) 



hup- 



2tt 2Q - 2p'_ 



n(p'_)rl(Q)tt{p'), 



(25) 

where vf- now denotes the Fermi velocity corresponding 
to the minority spin Fermi momentum kp-- A logarith- 
mic divergence appears as p' — > kp-- 

The above discussion describes the calculation method. 
We now need to calculate all the first order terms in the 
interactions g\ a p and g 2a p- The diagrammatic represen- 
tation of is shown in Fig. [21 The horizontal lines 
represent the electron being scattered by the Hartree- 
Fock potential of the Fermi sea (of scattering states). 
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Consider, for instance, the upper left diagram: an elec- 
tron, initially in state c& \ close to the Fermi level, passes 
through the barrier as a right-mover (a particle) and then 
interacts with the Fermi sea (on the positive z semi- axis). 
The electron is reflected (from a to b particle) while ex- 
changing momentum q with the Fermi sea. Finally, it is 
reflected by the barrier again, becoming a spin-up right 
mover with momentum p. A logarithmic divergence oc- 
curs if the polarized Fermi sea can provide exactly the 
momentum that is required to keep the electron always 
near the Fermi level during the intermediate virtual steps. 
Concerning the spin dependence of the interaction pa- 
rameters, we distinguish between g^, g 2 ^, which describe 
interaction between spin majority particles (that is spin- 
| on the right and spin- J. on the left of the barrier) and 
<7ii, (72J, j which describe interaction between spin minor- 
ity particles (that is spin-j on the right and spin-f on 
the left of the barrier). We use g±±, g 2 ± to denote inter- 
action between particles with opposite spins. According 
to the physical interpretation of the Feynman diagrams 
just given above, we always know on which side of the 
barrier the interaction with the Fermi sea (closed loop in 
the diagram) is taking place. 



p't Y q 01 




>'t k+q a k i a Pf 



Fermi level velocities v F ± for majority or minority spin 
particles, we write the diverging contributions to 6t^(p') 
as 



K, 



K F - 



K, 



rJO 


(oo I — 




AhvF— 




O — v' 




| 








- v' 


dQ 


(<72t — 


gu)r'(p')rUQ)t'Jp') 


4:hvp+ 




Q-p' 


dQ 




■9il)t\(p'-)rl{Q)r\(p') 


AhvF- 




Q-p'- 


dQ 


92±r>* 


{Q)W)r\{p') 


2hvF+ 


Q + 


Q--p'-p'- 


dQ 




(Q)t[{p'_)n(p') 


2hv F - 


Q + 


Q--p'-p'- 


dQ 


.92±f T 


(Q)r\(p%ip') 



2hv F+ Q + Q- -p' -p'_ 
Kf - dQ 52±rf(Q) n (p'_)4(p') 



2hv F _ Q + Q+-p> - p'_ 



(26) 



where Q± is related to Q as in Eq. (0J. In order to apply 
the poor man's renormalization method, it is preferable 
to transform the momentum integrations in Eq. i|26|) into 
energy integrals. In order to do this, we linearize the 
spectrum near the Fermi level as: 



hv F+ (Q - K F ^ 
hv F+ [p' - K F ^ 



nvF-(Q- - K F _) = e, (27) 
hv F -(p'_-Kp_)=e', (28) 



k 2« 



P't k+q a 



k_a 

I 2 



«2 



P't k+q a kjp pf 



FIG. 2: Feynmam diagrams for the first order contribution 
C/™ to the propagator 11911 . The scattering state is repre- 
sented by a double line, the a (b) particle is represented by 
a continuous (dashed) line. The loop represents the Hartree- 
Fock potential of the Fermi sea. The scattered electron ex- 
changes momentum q with the Fermi sea. 

It can be seen that the gi± terms are proportional to 
log \k F+ — k F -\, so they do not diverge. The logarithmic 
divergence would be restored in a spin degenerate system 
(/sf+ = kp— ). This can be understood from the diagrams 
in Fig. as follows: the electron with spin a is reflected 
by a polarized Fermi sea with spin —a. The momentum 
provided by the Fermi sea is 2k F - a , while the momentum 
required by the electron is 2k Fa . The g 2 ± terms produce 
logarithmic divergences that would not exist in the ab- 
sence of spin-flip scattering (i' = r' = 0). Introducing the 



where energy of the scattered electron is e' and the en- 
ergies e(e ) < are measured with respect to the Fermi 
level. The linearization is assumed to be valid within an 
energy range D around the Fermi level. The Q— integrals 
appearing in equation l|26|> can now be written as: 



dQ 



Q-P- 
dQ 

Q-p' 



-D 




— D 



de 



de 



K 



FT 



dQ 



Vp± 



de 

Q + Q±—p' — P- v F+ + vp- J_ D e — e' 



The scattering amplitudes with f (j) spin index are al- 
ways associated with the momentum p' (p'_). There- 
fore, we shall henceforth omit the momentum argument 
p' (p'_) of the scattering amplitudes. The divergent per- 
turbative correction, <$ij, is proportional to log(|e'|/-D), 



St i 



log^ 



im -flu) 

4hv F _ 



(ff2t - git) 

Ahv F + 



r T r T £ T 



(g2t - git) v ,/ {m -gu) * 

— ~ : ! ! — ~ l T I 



4hv 



F+ 



4hv F - 



6 



+ 



g2_L 



2h(v F + + vf-) 



(riV t t T + r'{ n t\ 



+ r' T *r T t' x + rf r' T t T ) . 



(29) 



For the calculation of St'^(p') and Sti(p'_), the prop- 
agators we need to consider are —(T T a p ,i(T)Cp, T ) and 
~(T T a p .[{T)c)pi respectively. The perturbation theory 
is analogous to that described above. In order to obtain 
t\{p'-) we consider the propagator —(T T a p ^(T)Cp, 

The logarithmically divergent perturbative terms can 
be dealt with using a renormalization procedure: we re- 
duce the bandwidth D by removing states in a narrow 
strip 5D near the band edge and work again the pertur- 
bation theory in the new bandwidth D — SD. The effect 
of removal of the band edge states must be compensated 
by adopting a new value of t-\ for the new bandwidth. 
Applying this procedure step by step yields successive 
renormalizations of t j . A differential equation is obtained 
by noting that the perturbation theory result, if + 5if, 
remains invariant as D is reduced: 

d6tf , 
*T + -q£ dD = 0. 

We introduce now a variable £ = log(D/D ) which will 
be integrated from to \og(\e'\/ Dq), corresponding to the 
fact that the bandwidth is progressively reduced from Do 
to \e'\ (which will eventually be taken as temperature: 
|e'| = T) and the scaling differential equations for the 
transmission amplitudes become: 



dt\ = (.921 ~ 9ij) 
ci£ Ahv F - 

(ff2T - git) 



4hv F+ 

92± 



2h(v F+ + 

+ r'fr^^+r'*^ ] , (30) 



dt' 



(.9: 



'21 



d£ 2hv F _ 

(ff2T ~ git) 
2hv F+ 

92± 



+ 



f r T*T 



+ 



h(v F+ + vf-) 



M'rrt + r\*r\t\ } , (31) 



d£ 4hv F + 1 1 111 

(921 - .9u ) r „* 



4Hvf- 

92± 



2h(v F+ +v F _) L'T i i 
+ rfri^ + rfrft ] , 



(32) 



d< i _ (.92 T -5it) l v , (g 2 ;-gu) . , 
t t J- + o>,„._ r i r i t i 



<K 2hv F+ ' ' 

.92 ± 



+ 



h(v F+ + v F -) 



2hv F - 

[r\*rih + r'*r\t[] . (33) 



In order to obtain the perturbative correction to the 
reflection amplitude r-\[p') we consider the propagator: 



G(r) = -(T r S PiT (r)4, ;T ) 



— e(p', |) p + p' + iO 



(34) 



In this case, there is a process where the incoming elec- 
tron from the left is reflected back by the Hartree poten- 
tial without even crossing the domain wall. The corre- 
sponding term comes from the Wick pairing term 

(bpA(r)bi^(r')) (4 iiT (r')S fel _ g , T (r')) (a k2+qA (r% A ) 

and gives a contribution to 5r^ (p') equal to 

.92T -giT / /\ i W\ 

The differential equation for r^(p') is 



+ 



+ 



.921 - gll r * 



Ahv F - 
92] - git 
4hv F+ 

92± 



2h(v F+ + v F -) 
+ 4*^4 + r{* t]t\ ] , 

and the differential equation for r'^(p') is 



(35) 



dr\ 



T . . g2T _ giT r * 



d£ 4hvF+ 



+ 



+ 



.921 ~ g u r„* 



[ r I*T*l + r l r i r T 



g2± 



2h{v F+ + wf-) 



L r T r i r T + r l r T r T 



92± 



2h(v F+ + v F -) 



1-(36) 



IV. FIXED POINTS 

The parameters of the model, which enter the scaling 
equations are: 



.92 t ~ git 
Ahv F+ 

921 - gu 
4hv F - 

92 ± 



2h{v F+ + v F _) 



= .9t - 
= .91, 
= 9± 



(37) 
(38) 
(39) 
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and the ratio vf-/vf+- The results can be presented in 
terms of the transmission coefficients, defined by 



T T = v -^-\u\ 2 

v F+ 1 11 ' ' | - ,u f \j logcrj 

T/ = |t' T | 2 , (41) | -is ^ — iog(r up ) 

2 



and the reflection coefficients 

= N 2 , 
= M 2 , 

K\ = ^|/ T | 2 

By definition, these coefficients refer to the respective 
currents divided by the incident current. 
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A. Insulator fixed points 




1 2 3 4 5 

| , 1 , 1 . 1 , 1 




-20 1 ' 1 ' 1 ' 1 ' 1 ' 1 

1 2 3 4 5 

log(D /| E |) 



FIG. 3: Logarithm of transmission/reflection coefficients ver- 
sus log^a. We may identify temperature T with |e|. The 
interaction parameters are: g-\ = gj, = 1, g± = 1.3 and the 
noninteracting domain wall model parameters are V — 0, 
vf-/vf+ = 0.8, vf+ = 1, A = 0.2. The dips are due to 
the sign reversal of the (small) scattering amplitudes. The 
long linear tails are analytically described in the text. At low 
temperature the system becomes a 100% spin flip reflector. 

We have made a numerical study of the scaling equa- 
tions. The non- interacting domain wall described in Sec- 
tion [n] provides the initial scattering parameters for our 
numerical scaling. Below we describe analytically the 
scaling behavior close to the fixed points we have found. 

For repulsive interactions (g-\ , ffj. , <7_l > 0) the system 
flows to insulator fixed points. For a moderate to large 
\/vp+ (larger than about 0.1) all the transmission am- 
plitudes, t a and t' a , vanish faster than any reflection am- 
plitude as T — > 0. We may then rewrite the scaling equa- 
tions neglecting the small transmission amplitudes. The 



FIG. 4: Logarithm of the transmission/reflection coefficients 
versus log(Z?o/|e|). We may identify temperature T with |e|. 
The parameters are:V = 0, (vf-/vf+) = 0.8, g-\ — gy = g± = 
1, vf+ = 1, A = 0.2. The reflection coefficients are all finite 
as T -> 0. 

scaling equation for r-f, for instance, becomes 

^ = 5T ( H 2 - 1 ) r T + 9l rlr{r\ + 2 9l _ ^|r' T | 2 r t , 

(46) 

where we used Eq. (|A10|) . The Wronskian relation ilA6(l . 
allowing for complex reflection amplitudes, shows that 

rJri+r T rJ*=0. (47) 

The charge conservation condition is satisfied solely by 
the reflections, 

1 = N 2 + — K\ 2 = h| 2 + — Kl 2 , (48) 

V F + V F + 

from which we easily conclude that \r^\ — \r jj at the fixed 
point. Then, Eq. I|46|l may be rewritten as 

-JF = ~ — ( 2 9± ~9t -91) r T r T 

«4 Vf+ 

= ^- (2. 9i -.g T -. gi )(l-|r T | 2 )r T . (49) 

VF+ 

Consider now the scaling equation J2fJJ). For in case 
of negligible transmissions we have 

dr' 

+ ff±(»fr;r T + ^1^1 ~ r \) ■ ( 50 ) 

The Wronskian relation l|A6(l . allowing for complex re- 
flection amplitudes, tell us that 

rtf + —r* r r\ = 0, (51) 
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and we recast Eq. i|50[) as 

dri 

-j7- = {gt + 9i-2g± 



VF- I 



J 1 2 



(52) 



In the derivation of 149(1 and (|52|l the only assumption 
made was that the transmission amplitudes are negligi- 
bly small. The reflection amplitudes may be, in general, 
complex and are still renormalized after the transmissions 
became negligible. 

Now we see that Eqs. (|49|) and (|52|l predict that the 
phases of the complex numbers rj , ri are unchanged dur- 
ing scaling. The two fixed points we may consider corre- 
spond to j-f approaching 0, or |rj| approaching 1 along a 
constant phase line in the complex plane. 

The situation |rj| — > requires 2g± — <7f — <7| > 
and, by charge conservation we have \r'^\ — > yJvp + /vF-. 
Upon integrating l|52|) with £ ranging from to 
log(T/Do)i the amplitude ri will vary from its initial 



Using the definition (|45|) for the re- 
flection coefficient, we write 



value rl to r'^ (T) 



i-n' , 



2(ST+Si- 2 9i) 



1 



TV, 



1-K' 



2(s T +9i-2SJ_) 



(53) 



If 2gj_ - g r - g l > 0, then K\(T) -> 1 as T -> 0. The 
domain wall becomes insulating. It reflects all incident 
electrons while reversing their spin. Therefore, such a 
DW may be considered as a perfect spin-flip reflector at 
zero temperature. In order to find the low T behavior of 
transmissions we put rp = rj_ = in equations H30fl - 133|) 
and obtain 



|t T |~|t' T |~|t' ; |~T 2 ^ 



(54) 



Figure [3] shows numerical solutions to the scaling equa- 
tions, where the system is flowing to this fixed point. 

In the regime where g^+gi — 2g± > we have IZ'-^ (T) — > 
0, TZ^(T) — > 1. So, the domain wall reflects all incident 
electrons while preserving their spin. From Eqs. I|3U|I - H33|) 
for the transmission amplitudes we obtain: 



I*tI 
141 
1*11 



y2gt 

J"2ffj. 



(55) 



For smaller values of \/vp+ (smaller than about 0.1) 
in the Hamiltonian (JIJ, the system flows to a fixed 
point, where rfj. vanishes faster than the transmissions 
and \r a \ — > 1. The transmission amplitudes still scale to 
zero as in Eqs. ||SSJ|. The scaling equation (|3*H|) for can 
be linearized in r'^ by neglecting the second order terms 
in t, t' , and considering that \r a \ — > 1, 



dri 



T' 



from which we obtain 



(56) 



(57) 



We see that the exponent for |r^(T)| is not bigger than 
the exponents in i|55|l . Now, in the scaling equation l|35(l 
for r-f we cannot neglect the terms containing transmis- 
sion amplitudes on the right hand side. One can eas- 
ily see that the g-\ term becomes 3|(|r-f| 2 — l)r-j-, which 
is of the same order of magnitude as the other terms. 
Consequently, the scaling behavior derived in Eqs. I|49|) 
and l(5^|) does not apply here, since it was assumed there 
that the transmissions were smaller than r' a . The be- 
havior of v a as T — ► can be found from the charge 

so 



conservation condition: = 1 — TZ\ — % — % 

1 — 1Zj ~ 7 1 min{2(g T +g i ),4g T } 



T _ -'T 

Such a situation is shown 



o 
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- i°g(R up ) 




2 4 6 

log(D /|E|) 



2 4 6 8 10 

log(D /|e|) 



If 91 + 91 - %9± = then both TZ\ (T) and !Z-\ (T) tend 
to finite values. Such a regime is illustrated in Fig.0] In 
this case, Eqs. (|30|) . I|31|) and i|33|l with constant reflec- 
tion amplitudes become a linear (in t^, t'p t'<) algebraic 
3x3 system. The eigenvalues of the matrix give three 
temperature exponents and each transmission amplitude 
will be a linear combination of the three powers of T. For 
decreasing temperature, there may be a crossover from 
one exponent to the other and the lowest exponent dom- 
inates as T — * 0. 



FIG. 5: Logarithm of transmission/reflection coefficients ver- 
sus log(Do/|e|). Wc may identify temperature T with |e|. Pa- 
rameters are: g-\ = 1, g± = 0.9, gi_ — 1.3 and vf-/vf+ = 0.8. 
For the initial noninteracting scattering amplitudes we used 
V = 0, v F + = 1, A = 0.05. 

in Fig. where t\ is seen to initially flow very fast to 
zero. The explanation is the following: for small A in 
Eqs. Q and ©, the noninteracting domain wall has 
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tf > 1, r-f > and tj < 1, rj. < 0. Also t' = r' is 
small. The scaling equation for t'< becomes 

dt' 

= (2 ffT r T + 2 9l _\ n \) r'fy + (2 5T r 2 - 2g x \r\r f l \) t\. 

The first term on the right-hand side is positive and much 
larger than the second one, so t\ tends fast to zero and 
disappears from the equations. The equation for tl is 



dti 



r',r[\)t' v 



The first term on the right is negative while the second 
is smaller because of small initial if. Then, initially 
grows as can be seen in Fig. [5] 



B. Transparent barrier fixed points 



i 

o 
o 

§ -1 

I -2 

h— "3 

-4 


_. -2 

55 
o 

3 _4 

z 

o 

I -6 
E -8 
-10 



l09(T up ) 

loaCTJ 



) 


0.4 


0.8 


1.2 1.6 










- iog(R dn ) 










iog(R up ) 










--- iog(R'J 















0.4 



0.8 1.2 

iog(D /|E|) 



Zero temperature fixed points corresponding to a 
transparent domain wall can be achieved when the in- 
teraction constants are all negative, i.e., for attractive 
electron interaction. Although we do not expect such a 
situation to occur in realistic physical systems, we de- 
scribe below the fixed points for the case V = in 
the model Hamiltonian (JIJ. For moderate to strong 
X/vp + in the model 0-©, the zero temperature val- 
ues 1 > = > |tf| depend on the initial param- 
eters. Smaller X/vf+ enhances if relative to t' a . The 
reflection coefficients vanish under scaling as powers of 
temperature. The corresponding exponents can be ob- 
tained after linearizing (for small reflections) the scaling 
equations (|3fci|) - The resulting 3 by 3 matrix contains 
the finite limiting values of the transmission amplitudes 
and its eigenvalues give the temperature exponents for 
the vanishing reflection amplitudes. Figure shows an 
example of this behavior. 

If some of the interaction constants are positive and 
the others negative, the situation becomes more complex. 
Below we describe several possible situations. 

1. The case <?f , g± > 0, g± < 

The system flows to the fixed point ff = rj. = — 1 
with all other amplitudes vanishing. The low-T behavior 
of the transmission can be easily found by inserting the 
fixed point reflections into Eqs. (|30|) - (|33|l : 

|t T | ~ T 9 T+9i , |t' T | ~ T 2ff T , \t[ \ ~ T 2g i . (58) 

The scaling equation for r+, neglecting second order 
terms in the scattering amplitudes, takes the form: 

dr' 

^l = (g 1 + gi -2g 1 _) r\ r\ ~ T^~ 2 ^ , (59) 

so that we must have g] + gi — 2g± > in order for 
r' ->0. 



FIG. 6: Logarithm of transmission/reflection coefficients ver- 
sus log(Do/|e|). Parameters are: <?f = —0.7, gy = —1.1, 
g± — —1 and vf-/vf+ = 0.8. For the initial noninteract- 
ing scattering amplitudes we used V = 0, vf+ = 1, A = 0.2. 
The transmission coefficients are all finite as T — > 0. 



2. The case g^ , gi < 0, g± > 

The system flows to the perfect spin-flip reflector fixed 
point |?*f| = ^Jvf+/vf- with all other amplitudes van- 
ishing, in accordance with the condition g^ +g[ — 2g±_ < 
derived earlier. 



3. The case g r > 0, g k < 

For a negative or small positive g ± , the system flows to 
a fixed point where \t'A = 1, rf = — 1. The wall transmits 
all spin-], particles with a spin-flip and reflects all spin-j 
particles. From Eqs. (|30|1 - (|33|) we see that the exponents 
for the transmission amplitudes are: 

|tf | ~ T 9r , \t\\~T 2 ^. (60) 

After linearizing Eq. I|36|) for in small amplitudes, we 
have 

dr' 

■^(ST-SxM r\~T^, (61) 

which requires g-\ > <?j_ for vanishing r^. If <?f — gj^ is 
small, the small quantities neglected in the right-hand- 
side of Eq. H6Ifl become important. Therefore, this fixed 
point holds for g^ — g± above some small quantity. Lin- 
earizing Eq. 1|34JI for rf in small amplitudes, we find 

^ ~ ~9l r l ~ 9iri n ~ T" 2 ^ , (62) 
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which tends to zero since gi < 0. For larger g± the 
system flows to the spin- flip reflector fixed point (\ri\ — > 

y/v F lJv F ~Z). 



4- The case < 0, gi > 



and down electrons, and the Fermi energy is E F 



h k F /(2m). Once the system becomes magnetized, the 
two new Fermi momenta, k F ±, satisfy the particle con- 
servation condition, 



fcp_ = 2k F => 



2, 



(67) 



The situation is analogous to the previous one. For 
negative or small positive g± the system flows to a fixed 
point where \t'^ \ = 1, — — 1 with all the others vanish- 
ing. The wall transmits all spin-f particles with a flip and 
reflects all spin- J. particles. From equations l|3T)jl - l|3*3'|) we 
see that the exponents for the transmission amplitudes 
are: 



(63) 



Linearizing equation (|36[1 for ri in small amplitudes, we 
have: 



drL t 



(64) 



which requires gi > g± in order for r'^ to vanish. If —g± 
is small, the small quantities neglected in the right-hand- 
side of Eq. 1)64(1 will become important. Therefore, this 
fixed point holds for g± — g± above some small quantity. 
For larger g± the sy stem flows to the spin-flip reflector 
fixed point (|rl — > v F+ /v F -). 



V. DISCUSSION AND SUMMARY 

Lateral ferromagnetic semiconductor wires with 
nanoconstrictions make it possible to achieve the limit 
of sharp domain walls^. It has been shown that the 
constriction itself does not cause significant reflection of 
the incident waves because it only produces a semiclas- 
sical potential^. We may estimate the parameter A of 
our model £p) by assuming that M(z) — Mo cos8(z)z + 
Mq s'm9(z)x with cos 9(z) = tanh(z/L)2&, where L is the 
width of the domain wall. We then find 



A 



J M n 



$m6(z)dz 



J Mo 

7T — L . 



implying that 



A 

VF+ 



JM 
(h 2 k 2 F+ /m) 



{Lk F+ ). 



(65) 



(66) 



The condition for the domain wall to be smaller than 
the Fermi wavelength is Lk F + < 2tt. The smaller Fermi 
wavelength is that of majority spin electrons, k F +- On 
the other hand, for small Lk F + the barrier becomes a 
poor spin-flip scatterer. The ratio v F -/v F + depends 
on the polarization degree of the electron system. Wc 
consider now a one-channel system. In a nonmagnetic 
system there is a single Fermi momentum k F for up 



and the spin-up and spin-down Fermi surfaces must cor- 
respond to the same energy, 



AE 



h 2 k 2 F _ 



AE 



(68) 



where AE/2 = J Mo is the Zeeman shift of the bands. 
From this we get 



k F ± AE 
k F AE F ' 

so that the ratio k F _/k F+ is 



VF- 
V F ^ 



1 - (AE/AE F ) 



k F+ v F+ 1 + (AE/4E F ) 
Inserting l|69|l into Eq. I|66|l we obtain 



A (AE/AE F ) 
vf+ [1 + {AE/AE F )Y 



(Lk+) . 



(69) 



(70) 



(71) 



The full polarization limit is k F - = and k F+ = 2k F , 
meaning that AE/AE F = 1, and then equation (17111 gives 



A 
v F+ 



0.79 Lk F+ . 



Typical values for a non fully polarized system are E F = 
90 meV and AE = 30 meV—. In this case we have 
v F ^/v F+ = 0.84 and equation l(7T|l gives 



A 

VF+ 



0.22 Lk F+ . 



Therefore, if Lk F + is smaller than about 2tt, the sys- 
tem can flow to any of the fixed points described above, 
especially the ones described in Section TlV Al 

The lateral quantization may produce several channels. 
The higher channels have larger Fermi wavelength and 
larger AE/AE F , so they can be in the spin-flip reflector 
fixed point. If a channel of high energy is fully spin polar- 
ized, then it corresponds to X/v F+ — 0.79 Lk F+ . But the 
possibility of inter-channel scattering arises. This could 
be due to the two following reasons: (i) electron-electron 
interactions (such would require a modification of our 
theory to allow for inter-channel scattering); (ii) the im- 
purity scattering. For the latter to be negligible we need 
the electron mean free path (not the transport mean free 
path) to be larger than the size of the constriction. 

In summary, we have studied the effect of electron- 
electron interactions on the transmission through a do- 
main wall in a ferromagnetic wire in the regime in which 
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the wall width is smaller than the Fermi wavelength. Ap- 
plying a renormalization technique to the logarithmically 
divergent perturbation, we obtained the scaling equations 
for the scattering amplitudes. The T = fixed points 
were identified corresponding to: (i) perfectly insulating 
wall (with or without complete spin reversal), and (ii) 
transparent wall. Both repulsive and attractive interac- 
tions were considered. We have estimated physical pa- 
rameters for a domain wall model which may be realized 
in physical systems. Such estimates suggest that realistic 
systems can display the behavior predicted in the vicinity 
of the fixed points we have found. 
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APPENDIX A: GENERALIZATION OF THE 
WRONSKIAN THEOREM TO SPINOR 
SCATTERING STATES 

The Wronskian theorem^ for (spin degenerate) scat- 
tering states in one-dimensional systems can be easily 
generalized to spinor states in a spin dependent scatter- 
ing potential. Let ipi(z) and ip2(z) represent two spinor 
scattering states with energies e\ and ti in the potential 
V(z). We assume V(z) to be a 2 x 2 real matrix, as is 
the case in the Hamiltonian Q , and consider a symmet- 
ric mass tensor rh (possibly position and spin dependent). 
Each spinor satisfies the Schrodinger equation 



d 1 dipi 
dz rh dz 

d 1 dip2 
dz rh dz 



ei-V 

€ 2 -V 



lp 2 



(Al) 
(A2) 



If we multiply the first equation (on the left) by tpK 2 ) 
the second equation by ip\{z) and subtract the two we 
obtain 



d_ 

dz 



t 1 # 2 lt 1 #1 



dz 



dz 



( £1 -62)^2, (A3) 



where the dot denotes the matrix (spinor) product (ip{ ■ 
ip2 = Ylcr rfiyO-ipZytr)- The expression in square brackets is 
a scalar function of z and would be proportional to the 
Wronskian of the functions ipi and tp^ in the case where 
the mass tensor reduces to a scalar. If the two states 



are degenerate (ei = £2), we conclude from i|A3(l that the 
expression in brackets is independent of the coordinate 
z: 

W(^,f 2 ) = • 1 • ^ -i>Uz) ■ I ■ y± = const. 

m dz m dz 

(A4) 

Since the potential matrix V is real, then ipl (or a l so 
satisfies the Schrodinger equation. 

The usefulness of the theorem expressed in equation 
(|A4|) is that it allows us to establish general relations 
between the scattering amplitudes, independently of the 
detailed form of the potential barrier. 

If we evaluate <|A4(1 for a pair of degenerate scattering 
states, say W(tp%. *, ipk-,i), the result must be the same 
for z < as for z > 0: 



which yields: 



z<0 



W(ipl v ip k -,l, 



z>0 



(A5) 



+ vrf(k)r' l (k-)+V-r\*(k)ri(k-)=0. (A6) 



Similarly, calculation of W(ip'l pip-k-,i) gives the rela- 
tion: 



«_ ^(fc)r i (A;_)+vt , f(fe)ri(fc_) 
+ v r|(fe)t i (fc_)+t!_ r\*(k)t\(k-) =0, 

and W(ipl p ip_f._ j i) gives the relation: 

w_ t]{k)r\{k) +v 4*(fc)r T (fc) 
+ v r^(k)t\{k)+v- r' T *(fc)t T (fc) = 0. 



A fourth relation can 



be obtained 



(A7) 



(A8) 
from 



i' i *(fc_)r i (fc_)+ U tJ(fc_)^(fc_) 
+ v rf(k-)t L (k-) + v- rj(fc_)tj(fe_) = 0. (A9) 

From W(if)k,t, V'fe-.l) we obtain the relation 

v r{(k-) =v- r\{k), (A10) 

and VF(^fe,T, ip-k-,t) gives 

»t i (L)=«-t t (t). (All) 

Considering a state and its conjugate, W(tpl M ,ipk,a) 
gives the conservation of the charge current, 

v = v_\m 2 + v\t\\ 2 + w|r T | 2 + v_\r\\ 2 , (A12) 

for <t =t, and 

v. = ^| 2 + v^\t[\ 2 + v-\ n \ 2 + v\r[\ 2 , (A13) 

for a =|. 



12 



APPENDIX B: FORMULATION FOR SPIN 
DEPENDENT ELECTRON EFFECTIVE MASSES 

Electron interactions such as (74 and <?2, which describe 
forward scattering between particles moving in the same 
direction, may produce renormalization of the electron's 
effective mass4£ The latter could depend on spin orienta- 
tion because the Fermi surfaces of spin-up and spin-down 
electrons are different. These effects can be taken into ac- 
count from the beginning by rewriting the Hamiltonian 
ijljl in a more general form, 

Ho = -#£^£ + nV6{z)-JM s {z)* M -nM{,)fi a , 

(Bl) 

where, in the kinetic energy term, we allow for a position 
and spin dependent effective mass tensor, rh{z). The 
tensor may take the form: 

*w=( m, o (2) ,4>)- (B2) 



with 



m-f(z) = TO + 0(— z) + m-Q(z) 



and 



m|(z) = m_0(— z) + to_|_0(z) , 



where &(z) denotes the Heaviside function. 

The appropriate mass values must be used in Eqs. (|2~|)- 
The expressions for the scattering eigenstates and 
transmission amplitudes given in the main text remain 
unchanged if we take into account that the velocities must 
be calculated considering the renormalized masses. 
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